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Abstract 

We give a gauge-invariant treatment of the angular momentum sum-rule for 
the proton in terms of matrix elements of three gauge-invariant, local composite 
operators. These matrix elements are decomposed into three independent form 
factors, one of which is the flavour singlet axial charge. The other two are interpreted 
as total quark and gluon angular momentum. We further show that the axial charge 
cancels out of the sum-rule. The general form of the renormalisation mixing of the 
three operators is written down and also determined to one loop from which the 
scale dependence and mixing of the form factors is derived. We relate these results 
to a previous parton model calculation by defining the parton model quantities in 
terms of the three form factors. We also mention how the form factors can be 
measured in experiments. 
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1 Introduction 



In the simplest parton model of the proton, its spin is carried by the spin of the valence 
quarks in such a way that their overall spin state is that of a spin-1/2 particle. In a QCD- 
improved parton model, however, such a picture must evolve with the renormalisation 
scale due to the splitting processes q — QCJ etc^ £it an arbitrary scale, some of the proton 
spin is then in principle carried by gluons and by strange quarks etc. Furthermore, it was 
realised some time ago by Ratcliffe |]J that the parton splitting processes generate orbital 
angular momentum. Thus, if the total spin of the proton is to have a scale-invariant 
meaning in the context of a parton model, then we must introduce quark and gluon 
orbital angular momentum components L q (x) and L g (x). Together with the usual quark 
and gluon intrinsic spin components Aq(x) and Ag(x), which are extracted from polarised 
Deep Inelastic Scattering (DIS), the sum of their first moments should be scale-invariant. 

The splitting functions for these four distributions were calculated by Hoodbhoy, Ji 
and Tang (HJT) to leading order. The form of the mixing of the first moments^ under 
scaling was found to be more general than originally envisaged by Ratcliffe 0. Impor- 
tantly, splitting functions do not depend on the choice of gauge used in their calculation. 
However, it is also claimed by HJT that the first moments of the angular momentum 
components, Aq, Ag, L q and L g , can be represented as forward matrix elements of four 
composite operators in proton states. Only one of these operators is gauge-invariant. 
This is the source of many difficulties, some of which are reported in a later paper |5|]. 
As an alternative, if we would define the four components of the proton spin in terms 
of forward matrix elements of gauge-invariant composite operators, we have a different 
challenge — there exists only one gauge-invariant operator which can measure total gluon 
angular momentum (Ag and L g in the parton model). This means one has to relate four 
parton model quantities to just three operators. This is discussed in refs. [[| [7|. Moreover, 
as far as we know, there exists no operator definition of L q (x) and L g (x) which scale in 
the manner of HJT. 

In this paper, we study quantities which have a precise operator definition. We discuss 
the gauge-invariant decomposition of the angular momentum current into three operators, 
thus restricting ourselves to quantities which can be interpreted as first moments only. 
We show that the forward matrix elements of the three gauge-invariant operators can be 
decomposed into just three Lorentz-scalar form factors. One of these form factors is the 
proton's axial charge which is measured from the flavour singlet component of the first 
moment of g\ in polarised DIS. Remarkably, we will show that the axial charge cancels out 
of the angular momentum sum-rule. We believe that this point has not been emphasised 
before. Since only two form factors appear in the angular momentum sum rule, we must 
conclude that the proton spin is only meaningfully decomposed into two pieces. We will 
interpret the two form factors as total quark and gluon angular momentum. Further, we 

1 Higher moments were calculated in refs. ft H. 
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need not interpret the axial charge as a measure of intrinsic quark spin. Indeed, because 
of the axial anomaly, the axial charge of the proton is, in the chiral limit, equal to its 
topological charge, which does not have an obvious interpretation as quark spin since it 
is defined as the matrix element of a gluon operator. An independent treatment of the 
anomalous suppression of the first moment of g\ (the "proton spin" problem) can thus be 
given in terms of universal topological charge screening ||. 

We also derive the renormalisation mixing of the three gauge-invariant operators. We 
show how, in terms of operator vertices, the axial anomaly affects two of the operators 
(one is the axial current) so that their sum is free of the anomaly. This is consistent with 
the fact that the axial charge does not appear in the sum-rule in terms of form factors. 
We have performed a one-loop calculation of the mixing of the operators. 

After presenting our calculations, we then show how they are related to the parton 
model calculation of H JT. To do this, we shall need to make an identification between the 
three form factors including the axial charge and the four parton model quantities. This 
is necessarily a Lorentz frame dependent statement, but which reproduces their results. 

Recently, it was proposed [0, |1(| that the two new form factors could be measured in 
hard, exclusive processes YP Xp such as Deeply Virtual Compton Scattering (DVCS) 
where X = 7. A section reviewing this possibility is included. 

The paper is organised as follows: In section ^| we review the decomposition of the 
conserved angular momentum current (AMC) into spin and orbital pieces, explaining 
why the gluon piece cannot be decomposed in this way gauge-invariantly. We include 
the gauge-fixing terms from the Lagrangian. Our main results are contained in sections [3] 
and [|: In section |3] we decompose the forward matrix elements of the three operators in the 
gauge-invariant decomposition of the AMC into scalar form factors; section |] contains the 
results of a one-loop calculation of the divergent pieces of the three operators, from which 
are derived the evolution and mixing of the form factors. In section |5| we discuss how to 
relate the form factors to parton model quantities, and how the result of HJT is obtained. 
In section |6] we review how the two new form factors can be extracted from a class of hard, 
exclusive processes such as Deeply Virtual Compton Scattering (DVCS). Finally, section [7| 
reviews attempts made in the literature to identify angular momentum components of the 
proton to matrix elements of non-gauge-invariant operators, highlighting the problems 
with such treatments. 



2 Decomposition of Angular Momentum Current 

In this section, we summarise the construction of a conserved angular momentum current 
(AMC) M tJ,I/X from the QCD lagrangian. As with the stress-energy tensor T^ v , we have the 
freedom to redefine M^ vX by adding on divergence pieces such that M fluX is still conserved. 
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As is discussed in ref. 0, it is possible in this way to separate the quark part of M^ uX 
into spin and orbital angular momentum (OAM) pieces. However, it is not possible to do 
the same for the gluon term in a gauge-invariant fashion. The reason is that the gauge 
field's spin and space-time indices coincide. Below, we extend the discussion in ref. [[J to 
include extra terms coming from gauge-fixing terms in covariant gauges and show that 
the decomposition of the AMC will depend on the details of the gauge-fixing procedure, 
so that even the ghost fields carry OAM. We will conclude from this that it is not sensible 
to define observables except in terms of the gauge invariant decomposition, even though 
this means we are left with only three operators measuring observable quantities instead 
of four. 

We begin with the divergence- free, symmetric stress-energy tensor T^ v . From it we 
can construct M fluX : 

M ^x = x i» T M^ (!) 

where the square brackets denote antisymmetrisation.0 This is conserved ( y d fl M fluX = 0) 
provided that T^ v is symmetric. It is always possible to redefine M fluX by adding the 
divergence of a function X mvX : 

M'» vX = M^ uX + d a X» auX . (2) 
M' is also conserved provided X is antisymmetric under [i <-> a as well as under v <-> A.Q 
If we now consider the covariantly gauge-fixed, BRST-invariant QCD lagrangian, 

C £gi + -£gf> 
£ gi = ^i>V>-ji^iV, (3) 
£ gf = -id^cD^c - WBAp + |5 2 , 

the corresponding symmetric stress-energy tensor is 

rpflU rpflU _|_ rpfMV 

T% = ^h^iD v ^ + F^F a v -gTC^ (4) 
Tg = -id {fM cD" } c-d^BA u} - g» v Cg. 

Writing M^ x = x^T^, one can show that 

A = ^ A >7 a 7 5 V> + tfrfx^i D x] ij + x ^F x]a F a " - x [u g x]fl C gi 
1 



- -d a 
4 



x^e^^ipjp^ip +EOM. (5) 



2 Our convention is that aH^ = a^b v - a"b^, while = c^b" + a v b^ . 

3 Note that a transformation T = T^ v + d a X^ va does not correspond to a transformation of Af . 
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Defining M' so as to cancel the divergence and equation of motion (EOM) pieces, we have 

m'/ vx = or A + o 2 M v + o£ [ v ] - x v 1 ^, 

M'f x = a>T g A f ]M (6) 

with 

Of = WiD v i>, (7) 
Of = F^FJ. 

We have in this way split the quark piece into what appears to be a spin piece, the axial 
current Oi, and an OAM piece XO2, while xO^ is the total gluon angular momentum, 
both spin and OAM. All the terms in M' gi are gauge-invariant. Matrix elements of M' gi 
vanish since it is a BRST variation of another operator. Also, forward matrix elements of 
the term propertional to x^ u g x ^ vanish; this will be shown in the next section. It is this 
decomposition which we will use to derive the scaling of the components of a nucleon's 
spin: It ensures that all quantities are explicitly gauge- invariant. 

Now consider what we get when we try to split the gluon component O3 into spin and 
OAM parts. Again, one can show that 

2 

_ F^ u A x] - F» a x [v d x] A a 

- A»x [u d x] B - id»cx [v d x] c - ix [v d x] cD»c - x [u g x] ^C 

+ d a [F" a x [u A x] ] + EOM. (8) 

We can redefine M' — > M" again by dropping the divergence and EOM terms. The new 
decomposition coincides with that of the canonical angular momentum current: Defining 
the rotation generator J tj , 

J ij = J (d 3 x)^, (9) 
if we choose the integral to be over 3-dimensional space at a given time, we have 
Jg) = J d 3 xM 0ij 

= J d 3 x [^a ij ^j + i^x {i d^ - E [t A j] + E k x [l d j] A k 
- A°x [l d j] B -id cx [i d j] c-ix [l d j] cD°c], (10) 
with E l = F° l and o % i = |[7*,7 :? ']. This may be written 

^)=E/ d ^^)(^Wa, (11) 
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where <p a are the canonical fields of the Kugo-Ojima canonical formulation of non-Abelian 
gauge theory [O], Tr((j) a ) are their conjugate momenta, and 5 iJ are their variations under 
a spatial rotation.^ 



One might be tempted to take the first four terms in eq. (§) and identify their matrix 
elements with quark spin and OAM and gluon spin and OAM respectively. There are 
good reasons for not doing so: 

1 . The other terms in eq. @ involving the gauge-fixing fields c, c, B are no longer a 
BRST invariant combination; the sum of their matrix elements cannot be expected 
to vanish. This means that we would have to talk about ghost OAM etc. in order 
to have a basis of operators which is closed upon renormalisation. 

2. The mixing amongst quark and gluon spin and OAM will depend on the gauge 
parameter, as observed in ref . . The observables whose scaling we wish to derive, 
on the other hand, are of course gauge invariant. 



3 Role of the Axial Charge 

We now develop further the gauge-invariant decomposition of the angular momentum 
current (AMC) in eq. @. First, we decompose the forward matrix elements of the oper- 
ators in the AMC into three Lorentz scalar form factors. One of these, the axial charge, 
appears twice with opposite sign, and so cancels out of the angular momentum sum-rule. 
This allows us to conclude that firstly the axial charge does not form part of the angular 
momentum sum rule and secondly that we do not need to interpret the axial charge as 
a measure of parton spin. The other two form factors have been discussed before || |7j]. 
however the relationship of these to the axial anomaly has not been emphasised. Fur- 
thermore, if the axial charge is to decouple from the sum rule in this way, we should be 
able to show this is reflected in the renormalisation of the AMC. To our knowledge, an 
explicit derivation of the mixing of the three gauge invariant operators appearing in the 
AMC has not yet appeared. Therefore, below we derive the general form of the mixing. 
In the next section we calculate the mixing explicitly to one loop. 

The axial charge clq is defined to be a Lorentz scalar form factor, or reduced matrix 
element: 

p } s) = Me^ Xp s p a Q} (12) 

We should be able to define all our observables in terms of such scalar functions, so let us 
now decompose forward matrix elements of all the remaining operators in the AMC (|B|) 

4 Each term in JiL satisfies the generator algebra for rotations. Note, however, that when the operators 
are renormalised, this is no longer true. Indeed, this even happens for the first term, the axial current, 
due to the axial anomaly. 
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except the divergence and the EOM, i.e. 0% x"\ O3 V' and g^ x x u ^C. To do this, we 
need to define the forward matrix element of an operator involving a factor of the co- 
ordinate x^. This is done in Appendix A. (See refs. ||, ||.) The result is that if B^ u {x) is 
a local composite operator, then the forward matrix element of x^B^°(x) is 



$ _1 (p) J d 3 x (p\x [i B l]0 ( 



x: 



a 



V 



5 J '1°(0) 



V 



(13) 



p'=p 



where $(p) is a wavepacket with no azimuthal dependence about the k-direction and 
where P = \{p + p') is held constant in the partial derivative w.r.t. A = p — p'. Before 
using this expression, we decompose the off-forward matrix elements of g^C, Of , Of 
and 



Of = Kp^i D v} i) 



(14) 



into scalar form factors. We write down only those terms which are linear in A and which 
are even in P; the latter constraint comes from symmetry under crossing, A — > A and 
P — > —P. The Lorentz structures must also have even parity: 



( P} s\g^C\p' } s) = C c (A 2 )M 2 g 



2 \W 



(p,s\0. 



p',s) = A q (A 2 )P^P-+^^P^e^iA a P pS(T + C q (A 2 )M 2 g^ 

+ ^p-PK^iAaPpSe + MD q (A 2 )e^iA aS/3 , (15) 



{p,s\Of(0)\p',s) = A g {A 2 )P^P^^^^P^e^iA a P pS<T + C g {A 2 )M 2 g 



{p,s\Of 



p',s) = A q (A 2 )P^ + ^^P^e v ^iA a P p s a + C q (A 2 )M 2 g 



(JLV 



Notice that the symmetric part of O4 is identical to that of Oi- Substituting these in the 
r.h.s. of we get 



d 



—1- 



dA u 



p,s 







—1 



—1 



—1 



OA,, 



d 
dAu 

d 
dA v 



p,s 



g» x £(0) 
Of(0) 



p',s 



p',s 



p=jy 



p=p> 



+ 



p,s 



p,s 



Of (0) 



p',s 



p',s 



p=p' 



p=p' 



0. 



2M P PP CT ' 

^P^°p p s a - MD q (p)<rW 8p , (16) 
^P^e^ Pp s a , 

2M P PP a - 
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There are further constraints on the five remaining form factors: We showed in deriving 
eq. (P) that 



0£ [ V ] = CT A + 0^x ui + divergence + EOM, (17) 
It is shown in Appendix A that the forward matrix element of the divergence term van- 



ishes. It follows from eqs. (|l4jlqJ17|) that 



B q (0) 
2D q (Q) 







a . 



Eqs. (0) and ( fl6|) therefore become 

(p, s or A (o) 

Of(0) 







dA,. 



d 

''dK 



0% x 



(0) 



p,s 



p,s 



V 



p,s 
«-> A) 



Maoe^sp, 



p=p' 



p=p' 



2M 



Sf3, 



(19) 



A) 



B g (0) 
2M 



The decomposition in eqs. (p!9|) is one of the main points of this paper: There can be 
no physically meaningful further decomposition of the proton spin in terms of forward 
matrix elements, since we only measure gauge-invariant quantities. Moreover, by adding 
the three equations together, one can see that the axial charge ao cancels from the sum- 
rule, so that only B q / g (0) appear in it. Thus, ao does not necessarily have to be interpreted 
as a measure of parton intrinsic spin. The best interpretation that one can give for B q / g (0) 
are the total angular momentum of the quarks and gluons J q / g respectively: 



B 9 (Q) 



J, 



(20) 



J, 



We refrain from decomposing J q / g further into spin plus orbital, since it is not clear that 
it is meaningful to do so. We shall discuss this further in section The sum-rule 



B q {0) + B g (0) = J q + J g 



(21) 



is assured by the non-renormalisation of the stress-energy tensor. This holds in spite 
of the presence of the axial anomaly because in eq. ([19]) it cancels out of the angular 
momentum sum- rule. 

In the next section, we calculate the scaling of the operators Oi^ 3 to one loop. Be- 
cause this is done in the language of operator vertices, as in the OPE analysis of DIS, it 
is important to show how the above conclusions translate into the language of mixing of 
composite operators. We shall argue that the axial anomaly affects not only the renor- 
malisation of 0\ but also of xC>2. To do this, we appeal firstly to some general properties 
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of local gauge-invariant operators, and secondly to the non-renormalisation of the AMC, 
which we have verified to one loop in section [|. We can thus write down the general form 
of the mixing of the operators of the AMC. 

First, we use the property that operators of the form xOi and O a when inserted in 
forward matrix elements mix with a block triangular matrix: 




7 -\ 





1. 1 



o b 



(22) 



This is assured because local, gauge-invariant composite operators with no factors of the 
co-ordinate x^ only mix with other such operators; the top right block of the matrix is 
therefore zero. Second, we note that because 0^ v is a symmetric tensor, it can only mix 
with the symmetric operators O3" and Of ■ 



it then follows that 0£ [ VJ 
.r" J and Of !yA + 0^ X x u \ when inserted in forward matrix elements. 



From eq. (FT 



can only mix with 0^ X x u ^ and Oj t!/A + 
Finally, given that the AMC is not renormalised, the columns of the mixing matrix should 
add up to one. With these contraints, the most general form of the matrix is as follows: 



/ 








< 1+X 










Z-X 1 + Z -Y 


V 




) 


B 


\ -Z -Z 1+Y 



0£ [ V 



(23) 



/ 



R 



where X = O(af), Y = 0(a s ) and Z = 0(a s ). The piece X is due to the anomaly, which 
appears at two loops. It cancels between 0\ and XO2 in the AMC. In terms of the form 
factors, this mixing matrix becomes 



a 

B q {0) 
B g (0) 




( «o 
B q (0) 



(24) 



In the above, we have taken the AMC current not to be renormalised. This we verify 
explicitly to one loop in the next section. 



4 One-Loop Evolution 

In this part, we present the mixing amongst the operators 0\, O2, O3 to leading order. 
First, we discuss the technicalities of operator mixing. The theory of mixing of gauge- 
invariant composite operators such as 0% is well understood. The extension to operators 
with a factor of x^ is straightforward. In particular, we show how operators of the form 
xOi can mix with operators of the form O a . 

Suppose we have a local composite operator B^ v (x) which mixes with a total derivative 
operator d a A flua (x). The forward matrix element of the total derivative operator would 
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vanish. However, an operator x x B flu (x) would correspondingly mix with x x d a A flua (x) 
which is not a total derivative, but which equals d a (x x A^ va (x)) — A tJ,uX (x). The first term 
has vanishing forward matrix element while the second is a local, composite operator with 
no factor x^ so gives an off-diagonal term in the mixing matrix. To see this explicitly, we 
represent the forward matrix element of x^ v B x ^{x) as (See Appendix A.) 



hm —i— — 

p'^p oA\ u 




I V 



(25) 



where is the momentum-space vertex associated with the operator B^, and the 
shaded blob is a four-point function, 1PI w.r.t. the p and p' legs. As an example, consider 
the UV divergence of a loop correction to T^ u : 

aT A] "(0 



d 



— i 




large r 



d 



dA w 




(26) 



The divergence piece takes the form, in dimensional regularisation, of 1/e pole pieces: 

F^(k, A, e) = F^(k, e) + iA a rr(k, e) + 0(A 2 ). (27) 

The 0(A) pieces, A a rf f correspond to mixing with a total derivative operator d a A^ va . 
The A differentiation acts on the blob, and then on the pole pieces. The first corresponds 
to mixing with another operator of the form xO^ while the second gives 

. d 



T x ^(k, A, e) = T x ^(k, e) - (i/ <-> A) + 0(A). 



(28) 



Thus, when considering the forward matrix elements, the mixing of an operator B^ v with 
a total derivative operator d a A^ va corresponds to mixing of B^ x x^ with A^ Xv \ 

Now we give the one-loop mixing of the operators in the gauge- invariant decomposition 
of the AMC. The basis of operators required is the following: 

Of" = V A > 7 «7 5 ^ 



Of = WiD v $, 
Of = F» a F a \ 



Of = -^iD u} ij. 
2 



(29) 
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These operators also mix with g^Os and g^Ei, where 



O, 



a/3, 



E x = itfjp^j, (30) 

but we shall not consider these since not only are they projected out in the definition of J 1 ^ 
in eq. but also the first derivatives of their matrix elements vanish in the forward limit. 
We do not require operators of lower dimension, since we are neglecting quark masses. 
Note that we require not just Oi^ 3 but also 4 , which is 2 with its indices symmetrised. 
The fact that 2 and O3 will mix with O4 does not spoil the non-renormalisation of the 
g.i. decomposition because of eq. (ITT 



Inserting the operators in the Green functions with two gluon or two quark legs, we 
calculate the pole pieces to one loop. Details are to be found in Appendix B. We find 
that the matrix required to renormalise the bare operators is 



(a 



B 



\ 



Of 

of 

or 



/ 



B 



ZiAo. 




V 



j)R, 

1 

: 








8 a s 1 
_ 3°^7 
8/~i a s 1 



1 9 A or A \ 



v 



or 

of 

or 



(31) 



R 



The angular momentum current is not renormalised to this order. Note also that the 
counterterms for O2 are the same as those for O4. Thus, the renormalisation of 2O2 and 
xO% follows that for the stress-energy tensor T^ v . At the two-loop level, we anticipate 
that this behaviour will be modified by the axial anomaly, where 0\ mixes with itself 
only jll : 

(Of A ) B = Z n (OfV (32) 
If the AM current is not to be renormalised, then xOi will have to mix with 0\. 



Using eq. fllTl) , when the operators are inserted in forward matrix elements, eq. 
becomes 

/ or A 



\ 



0&x x ^ 



( 







V 



8 a s 1 



1- fc^^J 

8ri a s 1 
3° F 4^7 





1 - ¥ft\ 



1 or A \ 

O^x^ 
{ Ofx x ^ ) 



(33) 



R 



If we now take the forward matrix elements of Oi, XO2 and xO% as defined by eq. 



we obtain, to 0(a s ), 



d 





dt ao = 




a s 






Jt Bq ~ 


47T 


d B 


a s 






dt B ° ~ 


47T 



~3 CfB * + 3 nfBa 



-C F B q 



; n f B 9 



(34) 
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5 The Parton Model Decomposition 



Since the angular momentum sum rule only involves two form factors and is independent 
of the axial charge and the anomaly associated with it, the best decomposition of the 
proton spin that we can have is 

Jq + Jg= \ (35) 

In the context of a parton model, we would like to be able to write 

^Aq + Ag + L q + L g = ^, (36) 

so that 

Jg = 7^Ag + L 9 , 
J g = Ag + L g . (37) 

This is not a natural thing to do since we would like to identify Ag with the axial charge, 
which multiplies a different Lorentz structure from J q / g in the form factor decomposition. 
In other words, such a decomposition is dependent on the Lorentz frame. 

Nevertheless, to make contact with the results of HJT |2j , we also relate the three form 
factors in terms of parton model quantities. As remarked above, the central problem is to 
relate four observables, Ag, Ag, L q and L g to just three form factors, not least because 
we can imagine each scaling in a distinct way. We will argue, however, that because 
it is possible to define Ag to be scale-independent, then we may determine the scaling 
of Ag, L q and L g gauge-invariantly in terms of the scaling of the three gauge-invariant 
operators in eq. (^). The problem is the same as that in polarised DIS, where there are 
two observables, Ag and Ag to relate to just one form factor, a . Below, we briefly review 
polarised DIS and the role of the axial anomaly, and how Ag can be constrained to be 
scale-independent in a certain class of renormalisation schemes. We then identify form 
factors with parton model quantities, and rederive the result of HJT || for the mixing 
of Ag, Ag, L q and L g . Thus we will have shown that an operator definition of these 
quantities does in fact lead to the HJT mixing. 

Let us begin with the first moment T\ of the polarised structure function g%{x, Q 2 ) as 
measured in polarised DIS. For large Q 2 , Ti factorises in the following way: 

r i ( e 2\ 
r 1 = j dx gi (x,Q 2 ) = ^[Cl iS (a s )Aq m + Cl(a s )Aq + 2n f C 1 g (a s )Ag\, 

Ag NS = -jrf^-lWft + Ag-J, (38) 
n f 

Ag = £(Ag, + Ag,), 

i=l 
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where (e 2 ) = -^jj^i^, Aq i: ^-li an d are the first moments of the polarised quark, 
antiquark and gluon parton distribution functions (PDFs) respectively, and C 1 (a s ) are 
the first moments of the co-efficient functions. Both the co-efficient functions and the 
PDFs are scheme dependent. At leading order, parton amplitude calculations show that 
the scaling of Aq and Ag are given by 




Os_ ( 
4tt { 3C F A) 




(39) 



where Cp = | and [3q = 11— |nj. The eigenvectors of this evolution are Aq, which does not 
evolve, and Aq—2nf^Ag. It then follows that there should exist a class of renormalisation 
schemes, the Adler-Bardeen (AB) schemes ||, where these two eigenvectors are the same 
to all orders. In such a scheme, it is the second of the eigenvectors which may then 
appear in the factorisation in eq. fl3"E|): By contrast, in the MS scheme, the gluon co- 
efficient Cg(a s ) vanishes to 0(a 2 ), so that the only flavour singlet piece in I\ is the 
contribution from Aq. The factorisation can also be expressed in terms of the operator 
product expansion (OPE). Since the only flavour singlet operator with the right dimension 
is the axial current, 

(p,s\or X \ P ,s) = Me^ xp sMt), (40) 

we have, for the MS scheme, a°(t) = Aq(t). That is, Ag does not contribute, while 
Aq(t) depends on the scale t because of the axial anomaly. However, in the class of AB 
schemes ||, the gluon co-efficient is C].(a s ) = — + 0(a 2 ), while C^as) = Cg 



a* 



— + 0(a 2 ), and we may identify |[| 



a (t) = Aq-2n f ^Ag(t). (41) 

Thus, the parton model interpretation of the axial current is that it measures both quark 
spin and gluon spin. The existence of schemes where this identification holds to all orders 
is assured by the Adler-Bardeen theorem fl3| . Aq is scale-independent to all orders, while 



the scaling of Ag(t) is determined by the scaling of the axial current. 

Now we make the parton model identifications according to the AB scheme interpre- 
tation: 

a = Aq -2n f ^-Ag, 

B q (0) = ^Aq + L q} (42) 
B g (0) = Ag + L g . 

The one-loop mixing derived in (|3~3| ) then gives us 

d \ 

j t A " = 



12 



dt qK ' 4tt 



d 
dt 



(L 9 (t) + Ag(t)) 



Aix 



8 1 2 
--C F (-Aq + L q ) + -n f (Ag + L g ) 

8 1 2 
-C F (-Aq + L q )--n f (Ag + L g ) 



(43) 



Thus dAg/dt remains undetermined. For this we require the two-loop evolution of 



15 



d_ 

di 



a°(t) = 7 a°(t), 



7 



-n f 



SO. 

2tt 2 



With the identification (f|l|), this implies 



d n/ 

-Ag(t) = ^[3C F A q + (3 Ag]. 
Combining fl4*3| ) and ([45]), we have the evolution equation, in agreement with ref. 



(44) 



(45) 



with 





I Ag ^ 




( 








^ 


( Aq \ 


d 


Ag 


_ a s 


3C F 


A) 








As 


~dl 


L q 


47T 


4 /-I 




— 


l n f 


L q 




V L g J 




V -\c F 


-li 




~3 n f J 


V ^ 9 / 



^-A q + Ag + L q + L g 



0. 



(46) 



(47) 



The general solution to these coupled equations is obtained in ref. 0, whence in the 
asymptotic limit t — > oo the partitioning of total quark and gluon angular momentum 
is \Aq + L q : Ag + L g = 3nf : 16. That result is a direct consequence of how the form 
factors B q / g (0) mix only with each other. This is the same partitioning between quarks 
and gluons of the first moments of the unpolarised PDFs. 



6 Measurement of Off-Forward Matrix Elements 



As was remarked above, of the three form-factors of operators appearing in the angular 
momentum sum-rule, only one, the axial charge, can be measured in polarised DIS. To 
measure the other two, one requires off-forward matrix elements of 0^ u and 0% u . As 
has been proposed in ref. J7|, [H|, these can be extracted from a class of hard, exclusive 
processes, generically 7*p — > Xp. For large Q 2 , these factorise into perturbatively cal- 
culable pieces, skewed PDFs and other non-perturbative pieces depending on the choice 
of the final-state species X. The skewed PDFs can be represented as off-forward proton 
matrix elements of a time-ordered bilocal operator product. The first moments of these 
are then equal to off-forward matrix elements of local composite operators. Factorisation 
theorems have been treated for Deeply Virtual Compton Scattering (DVCS) [T7], |TS 
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p p 

Figure 1: Definitions of four- momenta for hard, exclusive processes. 



7*p — > 7p, where the final-state photon is real, and also for meson production *y*p — > pV. 
For completeness, we review this work briefly, pointing out that we can measure proton 
spin components in the combinations B g (0) = J q and B g (0) = J g |7|, [TD[| . This supple- 



ments the isolation of Aq and Ag in polarised DIS and of Ag in open charm production. 



Consider then the process 7*p — > Xp, where we leave the species X unspecified. The 
amplitude /Ax for the process factorises as follows: 

M x (q,P,A) Q ^°J2Il [ duCi(Q 2 ,-,Z)fi/ P (u,Z,A*)D x/a . (48) 

The four-momenta are defined as in fig. |I|. The sum over i is over partons originating in 
the proton p; that over o is for, say, degrees of freedom going into the species X. The 
Lorentz-invariant kinematic variables are 



Or 



x 



-Q , 

Q 2 

2q-P 7 

q ■ A q2. 



2q ■ P 2P+ 

The unpolarised skewed PDFs, which we only need to consider, are given by 

f q/p (u,£,A) = |^e-^/ 2 ^-(p + iA|T(^)7Mr,0)^(0)) 
P+/ g/p K£,A) = / j^e- i( - u+ tW p+y ~ ^P + ^A \T(F +a (y-)V(y-, 0)F Q +(0)) 



(49) 
(50) 

(51) 



P A 

2 



P A 

2 



(52) 



where V(x, y) is the path-ordered string operator which renders a bilocal operator de- 
pending on two points x and y gauge- invariant: 



P(x, y) = P exp — ig / ds ■ A I . 



(53) 
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In the forward limit, A — > 0, these tend toward the usual PDFs with the correspondence 



/q/ P (w,f, A) -^2 /q/ p (u), 
f s / p (u,£, A) ^5 w / g/p ( w ). 



(54) 



The skewed PDFs have support for — 1 < u < 1, so their first moments w.r.t. u are 



du u /q/ p (u, £, A) 
1 

duf q/p (u,£,A) 



P + -A 

2 

P + -A 
2 



?/>7 + i £>+ ^(0) 



P-iA 
2 



F +a F a + (0) 



P--A). 



(55) 



Decomposing the matrix elements as in section |3| we can relate the form factors at zero 
momentum transfer B q (0) and B g (0) to the first derivatives of the first moments of the 
PDFs: 



where 



d f 1 



0A„ j-i 



duf g/p {u,£,A) 



A=0 



A=0 



M a /3 ' 



M 



(56) 
(57) 



B q (0) 
B 9 (0) 



(58) 
(59) 



Thus, if the spin vector s^-dependent part of the first derivative of the first moments of 
the skewed PDFs can be extracted from experiment, then this gives information on the 
total quark and total gluon angular momentum. 



7 Gauge-Non-Invariant Decompositions 

Here we discuss further the gauge-non-invariant decomposition of the AMC of eq. (Q), 
in particular the possible merits of relating this to the four components Ag, Ag, L q and 
L g in one particular choice of gauge and Lorentz frame. The natural choice for a parton 
model interpretation would be the axial gauge A + = and the infinite momentum frame, 
that is, the rest frame of the virtual photon in DIS in the limit Q 2 — > oo. Indeed, one 
finds attempts [|TI]] to relate moments of angular momentum PDFs, Aq(x), Ag(x), L q (x) 
and L g (x) to local, composite operators in the gauge A + = 0. However, one is henceforth 
restricted to a particular choice of gauge and quantisation procedure, viz. light-cone 
quantisation, which is inconvenient beyond low orders in perturbation theory. We also 
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discuss attempts to relate Ag to the forward matrix elements of the Chern-Simons current 
and review the problem which befalls this identification, namely the coupling of k^ 1 to 
unphysical pseudoscalar poles. 

To begin, suppose that we quantise QCD in the infinite momentum frame with the 
gauge A + = 0, i.e. use light-cone quantisation. In this picture, the only propagating 
degrees of freedom are the transverse gluon polarisations, A\, and we do not require any 
gauge-fixing fields or ghosts. Since the fields are quantised with equal x + commutation 
relations,^ then the canonical rotation generator is 



J* +) = J dx~d 2 x^ M +i \ 



(60) 



with 



Thus, 



1 



M i*x = -e^^yty + tyrfx^dflip - F^ U A X] - F» a x [u d x] A a . (61) 
2 



J (+) 



= J dx~ d 2 x ± [V2ipla ij ip + + iV2^ j + x [i d j] ^. 
- d + A [i A j] +d + A k ,x [i d j] A h 



(62) 



where ip+ = = -j^l This equation is of the same form as eq. (]TT ) . The 

difficulty here now is that the four operators are not even invariant under the residual 
gauge transformations which leave the condition A + = unchanged. In ref . |19 , however, 
definitions of Aq(x), Ag(x), L q (x) and L g (x) are constructed in the axial gauge which are 
invariant under residual gauge transformations. The first moments of these, J_ x Aq(x)dx 
etc., correspond to the four operators in eq. (|B"2"D, but with x^d^ replaced by x^V^, 
with being the derivative covariant under the residual gauge transformation. It is 
not altogether clear that those modified operators add up to give a conserved current. 
Moreover, one is restricted to a particular gauge, so if one encounters technical difficulties 
in perturbative calculations of the different contributions to the proton spin, as reported 
in ref. ||, one has no room for manoeuvre. 

The other feature of axial gauges is that the (+12) component of F^ U A X ^ is related to 
the Chern-Simons current k^. 



a 



k» = —e 
2tt 



s_ uvpo 



tr 



A u dpA a — A v A p A a 



(63) 



Indeed, putting A + = 0, we have k + = f^d + A [1 A 2 
and gauge-dependent statement. Some authors 

forward matrix element of k^ with — f^Ag in such a way that Ag is at least a Lorentz 



Note that this is a frame-dependent 
have tried to identify the 



12, 20, 21 



scalar form factor: From the anomalous Ward identity, which is not renormalised, 

d^-2n f Q~0, 



(64) 



5 Our convention is that x^ 1 = (x° ± x 3 ) 
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with Q = ^-ixF^F^ = and = ip^^ip, one deduces that — 2nfk^ is not 

renormalised, so that with 



(p,s\k^\p',s) = 2M/aJ(A 2 ) + s-AA%(A 2 ), 
(p,s\(jg-2n f k»)\p',s) = 2Ms»a° q (A 2 ) + s ■ AA%(A 2 ), (65) 



we have 

(p,s\^\p,s)=2Ms^a°(0) + 2n f a°(0)), (66) 



so that 



Aq = a° q (0) 



a. 



^Ag = oj(0). (67) 

This is true provided, as should be the case, that limA^o[2%Pc,(A 2 ) + p g (A 2 )} is finite. 

However, in covariant gauges, the presence of an unphysical massless pseudoscalar 
pole coupling to k^ means that as A — > 



2n f p g (A 



2\ 



1 

A 5 ' 

1 



p 9 (A 2 ) - -— . (6* 



In other words, the forward matrix element of the Chern-Simons current is singular, when 
treated non-perturbatively. In axial gauges, n- A = 0, this pole becomes a — pole, so that 
the Lorentz structure n M / (n • A) is present instead of A M / (A 2 ) in the limit A — > 0. This 
has been shown for the Schwinger model explicitly in ref. |23j and for QCD in ref. |p4}] . 



That said, for the gauge A + = 0, n + = 0, so no such pole is present in (p, s \ k + \ p, s); only 
for one choice of gauge and frame can the forward matrix element be well defined. 

In conclusion, defining observables in terms of matrix elements of gauge-non-invariant 
operators requires reference to a specific gauge, which is problematic. The matrix elements 
are also singular in the forward limit. Since the gauge- non- invariant approach is sufficient 
for our purposes, there is no need for further decomposition of the AMC. 



8 Discussion and Conclusions 

In summary, we have decomposed the conserved angular momentum current keeping it 
gauge invariant as far as possible. We have then decomposed the forward matrix elements 
of these operators into scalar form factors, and shown that the axial charge drops out of 
the angular momentum sum-rule. This is important because it provides evidence that 
the axial charge is not straightforwardly interpreted as a measure of parton intrinsic spin. 
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Indeed, the parton model interpretation which we have discussed in section |] is such that 
Aq and Ag appear in form factors multiplying two different Lorentz structures, one of 
which cancels out of the sum-rule. One should thus take the view that although it is 
meaningful to split angular momentum into quark and gluon pieces, a further split into 
spin and orbital pieces is frame dependent and therefore not a property of the proton 
itself. 

We have derived the mixing of the three gauge-invariant operators in the AMC, and 
calculated this explicitly to one loop. In general, we have shown how the axial anomaly 
cannot affect the renormalisation of the AMC, and this is reflected at the level of form 
factors by the decoupling of the axial charge from the angular momentum sum rule. 

Despite having argued that it is unnatural to do so, we have also derived the mixing 
of Ag, Ag, L q and L g by relating these to form factors of the gauge- invariant composite 
operators, and found results in agreement with what is obtained by taking the first mo- 
ments of splitting functions. The scaling of Ag has been determined from the scaling of 
the axial charge provided Aq is constrained to be scale-independent. This has enabled us 
to manage with just three operators. Thus, the mixing matrix of HJT can be derived in 
an explicitly gauge-invariant way in terms of local, composite operators. This means that 
we should not need to consider using gauge-non-invariant operators. 



Appendix A 



In this appendix we derive an expression for the forward matrix element of an operator 
containing a factor x M . This must be done with care otherwise we end up with a derivative 
of the momentum-conserving delta-function. To see this, consider the off-forward nucleon 
matrix element of an operator x^B^°: 



d 3 x (p 



x) 



V 



J d 3 xe l{p '- pyx x [l (p\B j]0 (0) 
i(2Tr) 3 -^- (5 3 (p-p')) (p 



P 



P 



(69) 



where {p \B(0)\p') is a Green function amputated w.r.t. its p and p' legs and where P is 
held constant in the derivative. We may fix this by replacing \p') with a wavepacket, 



d 3 A 
(2tt) 3 



$(p') \ p ') 



(70) 



such that Q(p') has no azimuthal dependence about the k-direction, perpendicular to the 
i- and j-directions. Then, integrating by parts, we have 



d 3 x (p 



x 



i B j]0 { 



d 



<9Ar 



v 



p 



p'=p 
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dAu 



V 



p 



(71) 



p'=p 



The first term on the r.h.s. vanishes since (p \B°i(0)\ p) = p°pp (p B p\ with (p B p) the 



reduced matrix element, so that we get a factor p^-J?—Q(p) which vanishes by construction 



i.e. the wavepacket carries no OAM. Thus we consider only the term 



$ 1 (p) / d 3 a; (p 



x 



d 



p 



p'=p 



(72) 



Notice that the l.h.s. reduces to (p\B(0)\p) if we remove the factor x M and talk about 
the matrix element of an ordinary, local, composite operator. 

Next, we show that matrix elements of operators such as d a (x l3 B(x)) vanish, so that 
the forward matrix element of the total derivative term in eq. @ can be neglected. We 
get 



d 3 x(p & y {x l3 B{x)) 



-(2 



7T 



d 3 A 

(2tt) 3 



(p - p') a 5 (p - p') 



d 



dA, 



i$(p'))(p\B(0)\ P ') 







+ Hp')-^(p\B(0)\p') 



(73) 



which is zero. 



Appendix B 

To evaluate the UV divergent pieces to one loop, we insert the operators Oi^4 into 
Green functions with either quark/anti-quark legs or two gluon legs, and calculate all the 
diagrams in figs. and|3|. The bare vertices, symbolised by ®, have the following Feynman 
rules, which also give the tree-level diagrams: 
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(76) 



Of: 





a. a 



or 





7,c /3, 6 a, a 
where P = i(p + p'), A = p — p'. 



2 1 



-p^p' v] g a p - p ■ p'g[ a gp } 
+P^g$p' p +p'^g; } p a 



-igf 



p^ gV}h gP]a _ g<*{^ gV}[l pfi] 
q{Vg"}l°lg'y]P — gP{l*gV}Wql\ 
fiV g v W g a ]l — g~l{V-g u W r a \ 



(77) 



(7f 



(79) 



We evaluate the - pole pieces of the Feynman diagrams and keep only terms of order 
A. For example, diagram Ela with ® = O2 gives us, in the Feynman gauge: 



C 



1 9 2 



e (4tt) 



C F - 



1 9' 



e (4tt) 



-o (? ?)t» + (1 \OZ u \ q) tree + - n 9» v (q l^ik) trce - (q \d a OT\ ?> t 



51) 



Adding up the diagrams in figs. ^] and |3| and their left-right mirror images, where appro- 
priate, we get, neglecting 0(e°), 0(A 2 ) and terms, 



or A 



1-Cp- 



1 g 2 



e (4tt) 



Of 



fiuX 



tree 



52) 
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(?ion?> 



B 



8 1 g 2 



+ 1 - Ci 



1 9 2 



e (4tt) 



<?iorig> 



tree 



(q\or\q) B 
{q\0TU) B 



8 1 g 2 



tree 



<<?|ori<?> 



B 



B 



finite 

2 
3 



-^f l T^(9\Or\9) tiee 



e (4tt) 2 



4 



(g\ong) 



3 3 



1 



e (4tt) 



Wis) 



tree 



53) 

54) 

55) 
56) 
57) 

58) 
59) 



Some of these poles come from the field-strength divergences of figs. 0b and |3|e-h and are 
absorbed by multiplying each quark and gluon external leg by Z^/ 2 and Z 1 / 2 respectively, 
with 



Z A 

Z^p 



1 



5 2 
3° F ~ 3*' 



1 9 2 



e (4tt) 



1 9 2 



e (4tt) 



(90) 



Writing 



(q\o b \q) r 
(9\0 B \g) R 



Z<p (q\0 B \q) B , 

ZA(g\0 B \g) B , 



we have 



(9 



01, B 

2 , B 



9>a 



(9|Oi, B |g) tree 

8 1 g 2 

(q \o 2 ,b\ q) tIce - 3^-7435 I°4,b| 9) 



tree 



1 9 

3 F e(47r) 2 



9) 



(9|04,s|g) tree 

{qpAAq) t 



1 s 2 



1 — —Cp 

\ 3 e(47r) 2 7 

finite 

h\ik? {9p * Ag) * 



(91) 

(92) 
(93) 

(94) 

(95) 
(96) 
(97) 
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3'"' e (4?r) 

(^|0 4)B |^ = (^|0 2 , B |^. (99) 



|0 3 ,B|p) fl = l--n /77 — ^ ^|0 3) B|^} tree (91 



The renormalisation matrix required to remove these remaining poles is then given by 
eq. ©. 
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Figure 2: One-loop diagrams 
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Figure 3: One-loop diagrams 
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